Goos-Hanchen induced vector eigenmodes in a dome cavity 
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We demonstrate numerically calculated electromagnetic eigenmodes of a 3D dome cavity resonator that owe 
their shape and character entirely to the Goos-Hanchen effect. The V-shaped modes, which have purely TE 
or TM polarization, are well described by a 2D billiard map with the Goos-Hanchen shift included. A phase 
space plot of this augmented billiard map reveals a saddle-node bifurcation; the stable periodic orbit that is 
created in the bifurcation corresponds to the numerically calculated cigenmodc, dictating the angle of its "V" . 
A transition from a fundamental Gaussian to a TM V mode has been observed as the cavity is lengthened to 
become nearly hemispherical. © 2008 Optical Society of America 
OCIS codes: 140.4780, 260.5430, 260.2110, 230.5750 



The Goos-Hanchen (GH) effect, 1 which typically oc- 
curs upon reflection from a dielectric interface or layered 
structure, is an apparent parallel translation of the re- 
flected beam away from its geometrically expected posi- 
tion. A simple treatment 2 shows that the size of the GH 
shift is proportional to the derivative of the reflection 
phase with respect to the angle of incidence. In addition 
to shifting beam position, the the GH effect can manifest 
itself in alterations of differential cross sections, 3 of laser 
mode dynamics, 4 and of mode spectra. 5 

Here we present and analyze a class of eigenmodes in a 
non-paraxial dome cavity, the "V modes" , showing that 
the modes are straightforwardly explained by a stable 
periodic orbit in a billiard map that is augmented by the 
Goos-Hanchen effect. A ray analysis of the cavity geom- 
etry without taking the Goos-Hanchen shift into account 
fails to predict these modes. To the authors' knowledge, 
this Letter presents the first demonstation of the GH 
effect, or any nonspecular effect, 6 creating entirely new 
eigenmode patterns. 

Our model dome cavity consists of a perfectly conduct- 
ing concave spherical mirror, Ml, facing a planar thin 
film (distributed Bragg reflector), denoted by M2. The 
electromagnetic eigenmodes are calculated by the basis 
expansion methods which are fully described in Ref. 7. 
We first described the V modes in Ref. 7, but no physical 
explanation for their appearance has been put forward 
until now. 

To provide quantitative results, we fix some parame- 
ters at the outset, cf. Fig. QJa): The radius of Ml is R = 
60 fim and the length of the resonator, measured from 
the front surface of M2 to the zenith of Ml, is L = R—Zi, 
where we will focus on the range 0.3 < z\ < 0.7 /j,m. 
This means the resonator is nearly hemispherical, and 
increasing mirror offset z\ makes the cavity more parax- 
ial. The layer structure of M2 is IAAAA{BA) 22 where 
/ denotes the vacuum region inside the cavity, each A 
denotes a layer of n = 3.52 material (AlGaAs) of quar- 
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Fig. 1. (a) Diagram of the simple ray model for the V 
mode, in a dome cavity, defining the polar angle 9 and 
the mirror labels Mi, M2. (b) Gray-scale plot of a VTEq 
mode at mirror offset z\ = 0.4 fim. Parts of the dome 
surface are outside the figure for space reasons, (c) An- 
gular spectrum of this mode. 



ter wave thickness at 570 nm, each B denotes a layer of 
n = 3.0 material (AlAs), also of quarter wave thickness 
at 570 nm, and O denotes the vacuum region outside the 
cavity. 

Axial symmetry allows one to label all cavity modes 
by a total angular momentum quantum number to; the 
complex electric and magnetic fields E p , E^, E z , B p , 
Bfi, and B z (in cylinder coordinates) are proportional to 
cxp(iTO</>) . In this paper, we present data only for modes 
with m = ±1. The vectorial eigenmodes typically found 
in an axially symmetric cavity are in general neither TE 
nor TM. Here, TE refers to fields whose angular spec- 
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trum contains only s-polarized plane waves with respect 
to the Bragg mirror; TM denotes a p-polarized plane 
wave decomposition. 

Fig. \V[b) shows a V mode having wavelength A = 
755.44 nm at z\ = 0.4 ^m. The contrast-enhanced plot 
shows the y-z slice of the physical E x field (z is the dome 
axis). Azimuthally, the electric field here is proportional 
to exp[i(±0 — tot)], which may be combined to make the 
"linearly" polarized versions cos/sin((/>) exp(— iut). The 

V mode is almost completely TE polarized. Fig. [TJc) 
shows the s-polarized (TE) polar angle distribution of 
plane waves, \ip s (9k)\ 2 , for this mode; 9k is the angle of 
incidence of a plane wave with respect to the cavity axis. 
The p-polarized admixtures are negligible on the scale of 
the plot. 

Averaging 9k against this distribution yields (9k) = 
53.6°, which one can identify as the "opening angle" of 
the V-shaped beam. Similar modes have been observed 
which arc entirely TM. We also found TE modes which 
exhibit a transverse node in each side of the V. This leads 
us to introduce the notation VTM™ and VTE™ for these 

V modes, where m is the total angular momentum and 
v is the number of nodes in the transverse wavefunction 
of each side of the V. The mode in Fig. Q] is then VTEj. 

We now turn to the explanation of these V modes in 
terms of the Goos-Hanchen effect: The GH shift along 
the reflecting surface, of an s- or p-polarized beam is 2 
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where k = 2tt/X and 4> s / p is the plane wave reflection 
phase of the surface. A simple self-retracing periodic or- 
bit which reflects perpendicularly from M2 can exist at 
9k = x > if the GH shift equals the geometrically 
determined segment, that is if 

(2) 



A s / P (x; k) = 2z! tanx, 



where the right hand side is seen to be the length of 
the horizontal segment of the orbit shown (bold line) in 
Fig. IHa). 
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Fig. 2. (a) Comparison of (9 k) with x as a function of z\ 
for a VTEj mode. The \s calculation takes its k values 
from the actual modes, (b) LHS, RHS of Eq. j2|) (using 
A s ) with respect to x (deg), with fixed A. The solutions 
corresponding to (s)table and (u)nstable orbits alternate. 

Since z\ is given physically and k may be coarsely 
chosen, Eq. |[2|) is in practice solved for Xs or Xp to pre- 
dict the angle of the V modes. The prediction for the 



VTEj mode shown in Fig. [T] is \ s = 54.2°, in excellent 
agreement with the average angle (Ok) = 53.6° obtained 
for the numerically calculated eigenmode. Numerical ex- 
perimentation shows that Xs/p increases to near 90° as 
z\ — > 0+. Fig. HJa) shows good agreement between Xs 
and (Ok) for a VTEj mode followed from z\ = 0.3 to 
z\ = 0.49. (Following the VTEj mode pattern over this 
range required "hopping" across mode anticrossings at 
z x = 0.46 and zi = 0.48.) 
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Fig. 3. Billiard map with GH shift for TE polarization 
at z\ = 0.55, showing creation of the stable V-ray island 
and separatrix. 

To characterize the conditions under which the V 
modes are stable, it is useful to investigate the phase- 
space of the classical ray dynamics. The ray orbit shown 
in Fig. QJa) is a periodic orbit (PO) of a billiard map 
which obeys specular reflection at the top mirror M2, 
but at Ml receives a GH shift, along the Ml surface, 
the length of which is given by Eq. (Q]). To specify this 
map, we follow a ray through multiple reflections, and 
just before the n-th bounce at mirror M2, note its polar 
coordinate 0^ as shown in Fig. [Ha), as well as the com- 
ponent of the (unit) momentum in the direction. 
The billiard map T maps this pair of ray properties from 
one reflection at M2 to the next one: 



r:(pW *W)~(p<^V B+1) ), 



(3) 



A plot of nine trajectories, each of 1000 bounces, of the 
map T for z\ = 0.55 and s-polarization (TE) is shown 
in Fig. [3l Such a phase space plot is very much like the 
Poincare section of a dynamical system and may be ana- 
lyzed as such. The island structure surrounding the sta- 
ble 9 = PO manifests itself in wave physics as the fam- 
ily of paraxial Gaussian modes. The separatrix indicates 
an unstable PO at 6 rs 37° and the second island indi- 
cates a stable PO at 9 ~ 49°. As shown in Fig.^b), these 
two POs are created at the same point (this is called a 
saddle-node bifurcation) as z\ is lowered to about 0.6, 
and move apart in 9 as z\ is further decreased. As z\ 
decreases, the 9 = island shrinks, while the new island 
grows. As it grows, one or more existing wave modes 
may become distorted until one of them becomes the V 
mode; the V modes live inside the new island. If the GH 
effect is turned off in the billiard map ([3]), there is no 
new island, and the only PO that exists is the 9 = 
orbit. This leads to the important conclusion that only 
a ray map augmented by the GH effect predicts the sta- 
ble cavity modes correctly. This analysis can be applied 
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equally to the s- or p-polarized V modes, by inserting 
the appropriate reflection phase (j)s/p m the map T. 

At sufficiently large mirror offset Z\, we should be 
able to make contact with Gaussian-beam theory. The 
resonator modes may be characterized by their trans- 
verse field patterns, which are described in Rcfs. 7—11. 
In Rcfs. 7, 10, 11, a dielectric mirror provides a form bire- 
fringence which causes the transverse mode patterns of 
many typical cigcnmodes to be symmetric and antisym- 
metric mixtures of specific pairs of vectorial Laguerre- 
Gauss modes, causing either TE or TM polarization 
to dominate the plane-wave spectrum. Starting from 
paraxial theory, these "mixed" modes can be denoted 
by the transverse order N , the angular momentum m, 
and whether or not they arc quasi-TE or quasi-TM: 



QTE(M)™ = LG 



(JV — |m- 
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for \m\ > 1, TV - |m+ 1 £ {0, 1, . . .}, and N-\m-l\ e 
{0, 1, . . .}. Here cr^ is the Jones vector for right/left cir- 
cular polarization and LG denotes the Laguerre-Gauss 
wave pattern. 10 A detailed discussion of the mixed modes 
is given in Rcf. 11; they, along with pure Gaussians oc- 
curing for |m| = N + 1, comprise the background of 
cigcnmodes against which the V modes appear. 
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Fig. 4. Angular spectra showing the creation of the 
VTMj mode from the fundamental Gaussian as z\ de- 
creases. Dashed: \ip s \ 2 , solid: \ip P \ 2 - 

We numerically observe that V modes evolve from 
some existing mode of this paraxial type as Z\ is re- 
duced and the stable island grows. As an example of 
this metamorphosis, a VTEj mode was created from 
the predominantly TE mixed mode having the tran- 
verse E field, in the paraxial limit, given by QTE4 = 
LG°,er + +LG'j ! o-~. On the other hand, as shown in Fig. 0] 
a VTMj mode was created from the fundamental Gaus- 
sian mode TEMoo, which originally has circular polariza- 
tion. The polarization-dependent stabilization provided 
by the GH effect thus overcomes the polarization state 
expected for the paraxial fundamental mode in an ax- 



ially symmetry cavity. This non-paraxial nature of the 
V modes underlines the predictive power of the Goos- 
Hanchen picture. 

Although the numerical results presented here were 
limited to |m| = 1, we have observed VTM§, VTEq, and 
VTEj, and expect V modes to generally exist for |m| > 
1. Other dielectric stack designs than the one described 
here have also been observed to induce V modes. 

We note that patterns reminiscent of the V modes ap- 
pear in an earlier work by Blumel et al. 12 which discusses 
scalar waves in a purely two-dimensional circle with a di- 
electric boundary across its diameter. This roughly cor- 
responds to the limit z\ — > in our model, and neither 
the question of the GH effect nor the z\ dependent open- 
ing angle of the V modes are addressed in that work. 
As pointed out in Ref. 12, semiclassical techniques that 
may be appropriate for smoothly varying media or sim- 
ple boundary conditions (e.g., Dirichlet) need to be ex- 
tended in the presence of dielectric interfaces. The results 
presented here were obtained in a simpler and physically 
appealing way, by adding nonspecular corrections to the 
classical ray mapping T. 

The V modes demonstrate that small corrections to 
specular reflection, such as the Goos-Hanchen shift, can 
have large qualitative effects. Billiard models augmented 
by nonspecular corrections possess new phase space 
structures which correspond to the new modes. The V 
mode in particular provides a potential lasing mode with 
two or four entrance/exit channels, yet it economically 
uses only one curved surface. Interestingly, the spread- 
ing angle of each limb of the V modes is smaller than 
predicted by a Gaussian-beam calculation based on the 
path length of the V including GH shift. In future work, 
we will investigate to what extent a stability analysis of 
the ray map T is able to explain this phenomenon. 

This work was supported by NSF Grant ECS-0239332. 
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